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ABSTRACT 

Context. Various studies have emphasized variations of the charge state and composition of the interstellar polycyclic aromatic hy- 
drocarbon (PAH) population in photodissociation regions (PDRs). These changes are expected to impact the energetics and chemistry 
in these regions calling for a quantitative description. 

Aims. We aim to model the spatial evolution of the charge and hydrogenation states of PAHs in PDRs. We focus on the specific case 
of the north-west (NW) PDR of NGC 7023, for which many observational constraints are available. We also discuss the case of the 
diffuse interstellar medium (ISM). 

Methods. The physical conditions in NGC 7023 NW are modelled using a state-of-the-art PDR code. We then use a new PAH 
chemical evolution model that includes recent experimental data on PAHs and describes multiphoton events. We consider a family of 
compact PAHs bearing up to 96 carbon atoms. 

Results. The calculated ionization ratio is in good agreement with the observed ratio in NGC 7023 NW. Within the PDR, PAHs 
evolve into three major populations. We find medium-sized PAHs (50 < Nc ^ 90) to be normally hydrogenated, while larger PAHs 
(At <; 90) can be superhydrogenated, and smaller species (At ^ 50) are fully dehydrogenated. In the more diffuse gas of the cavity, 
where the fullerene C60 was recently detected, all the studied PAHs are found to be quickly fully dehydrogenated. PAH chemical 
evolution exhibits a complex non-linear behaviour as a function of the UV radiation field because of multiphoton events. Steady state 
for hydrogenation is reached on timescales ranging from less than a year for small PAHs, up to 10 4 years for large PAHs at Av = l. 
Critical reactions that would need more studies are the recombination of cations with electrons, the reactivity of cations with H 2 and 
the reactivity of neutral PAHs with H. 

Conclusions. We developed a new model of PAH chemical evolution based on the most recent available molecular data. This model 
allows us to rationalize the observational constraints without any fitting parameter. PAHs smaller than 50 carbon atoms are not 
expected to survive in the NGC 7023 NW PDR. A similar conclusion is obtained for the diffuse ISM. Carbon clusters turn out to be 
end products of PAH photodissociation, and the evolution of these clusters needs to be investigated further to evaluate their impact on 
the chemical and physical evolution of PDRs. 

Key words, astrochemistry - ISM : molecules - molecular processes - photon-dominated region (PDR) - Methods: numerical 



1. Introduction 



<^> Polycyclic aromatic hydrocarbons (PAHs) are now widely ac- 
C*~) cepted as the carriers of the ubiquitous aromatic infrared bands 
1—1 (AIBs) observed at 3.3, 6.2, 7.7, 8.6 and 11.3 yt/m and this has 
motivated a lot of studies on the properties of these species (see 
|Joblin & Tielens 201 1 for a recent compilation of papers on this 
subject). Interstellar PAHs can contain up to 20% of the cosmic 
carbon and are the smallest dust particles by size. As such, they 
play a major role in the physics and chemistry of photodissoci- 
ation regions (PDR) through the absorption of UV light (Joblin 
et al.||1992) , the photoelectric heating ( [Bakes & Tielens||1994jl 



X 

C3 



and may also be involved in the formation of H2 (Habart et al. 
[2003] |. 

The energetics of PDRs is well understood in its main lines. 
However, it has been known for more than a decade now (see, 
e.g. |Hollenbach & Tielens|1999| l that the excitation temperatures 
observed for some molecular tracers are higher than predicted 
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by PDR models. This issue was raised again recently by several 
authors who compared the results of the most recent PDR mod- 
elling tools with observational data from the last generation of 
spatial facilities (Spitzer, Herschel). These authors reported diffi- 
culties to account for the observed emission in H2 and high-J CO 
rotational lines ( Joblin et al.|2013a, in prep.[|Habart et al.|201 1[ 
|Goicoechea et al.|201 \\ . |Habart et al.| ( |201 \\ concluded that this 
is due to our poor understanding of the energetics of PDRs. The 
main heating mechanism at the surface of PDRs comes from the 
thermalisation of hot electrons generated by the photoelectric ef- 
fect on PAHs and very small grains. The photoelectric efficiency 
depends on the grain charge and is expected to decrease when 
PAHs are positively ionized (Hollenbach & Tielens 1999). This 
theoretical prediction has been recently supported by the analy- 
sis of observations (Okad aet al.| submitted). 

The question of the exact evolution of PAH properties with 
physical conditions has been addressed recently by a combina- 
tion of observations and theoretical modelling. By use of the 
Infrared Space Observatory (ISO) and the Spitzer telescope, the 
spatial evolution of the mid-IR spectra was observed in several 



Article number, page 1 of 



16 



regions. Rapacioli et al.| ( |2005| ) and |Berne etaE] ( |2007| l used 
global decomposition methods to analyse the data and proposed 
an evolutionary scenario in which a population of very small 
grains evaporates under the action of UV photons, releasing neu- 
tral free PAHs, which are subsequently photoionized. In addi- 
tion, the recent detection of fullerene Ceo in reflection nebulae 
( |Sellgren~e t al. 2010|l may be interpreted as a new step in this 
evolutionary scenario (Berne & Tiele nsp012| l. 

In parallel with these observational studies, several models 
were developed to describe the chemical evolution of PAHs. 
Bakes & Tielens (1994) focused on the charge evolution in or- 
der to compute the contribution of PAHs to the gas heating by 
the photoelectric effect. A more comprehensive model was pro- 
posed by Allai n et aL] ( [T99 6a ), which includes the photodestruc- 
tion of PAHs exposed to UV photons, and provides insights into 
the stability of interstellar PAHs, as a function of their size. 
|Le Page et al.| ( |2001| l proposed a new model dedicated to the 
charge and hydrogenation states of PAHs in the diffuse interstel- 
lar medium (ISM). Based on these models, Vis ser et aL] ([2007 ) 
built a comprehensive model of PAH evolution in protoplane- 
tary disks, which includes the charge and hydrogenation states. 
The authors also described the photodestruction of the skeleton 
of PAHs including multiphoton events that are frequent in these 
high-UV irradiation conditions. 

In this paper, we present a new model that is inspired from 
; Page et al. (2001 1, and includes multiphoton events for all 
photodissociation processes. This model takes advantage of the 
latest available molecular data from both theoretical calcula- 



Table 1. Input parameters of the Meudon PDR code ( [Le Bourlot et al. 
|2012) used to model the physical and chemical structure of NGC 7023 
NW PDR. 



tions ( Malloci et al.|2 007 , 2008) and experimental measurements 
(|Joblin et a l. 2013b, in prep. Biennier et al.|[2006{ |Betts et al.| 



2006 ). We apply our model to the prototypical north-west PDR 
(hereafter NW PDR) of NGC 7023 in an attempt to (i) identify 
key processes in the evolution of the charge and hydrogenation 
states of interstellar PAHs in order to provide guidelines for fun- 
damental studies and (ii) provide a quantitative description of the 
evolution of the PAH population, which is necessary for mod- 
elling the role of PAHs in the physics and chemistry of PDRs. 
Furthermore, the obtained results could guide the spectroscopic 
identification of individual PAHs for instance by matching some 



of the diffuse interstellar bands (see, for example, |Salama et al. 
|2011| l. This motivates additional calculations on the hydrogena- 
tion and charge states of PAHs using a grid of conditions that 
include typical conditions for the diffuse ISM. 

This work is organised as follows. The studied PDR and its 
physical conditions are presented in Sect. [2] In Sect. [3] molec- 
ular processes relevant to this study are reviewed and discussed. 
Our numerical model is presented in Sect. [4] and its results in 
Sect. [5] The sensitivity of the results, both to molecular data 
and astrophysical conditions, are investigated in Sect. [6] Conse- 
quences of our results on the direct detection of individual PAH 
species and on the fate of carbon clusters are discussed in Sect. [7] 



2. The astrophysical environment 

Our study focuses on the reflection nebula NGC 7023 that has 
been widely studied at many wavelengths (e.g. |Rogers et al 
T995]|Fuente et al |19961|Gerin et al.|f if9"8"t|Fuent e et al.||1999[ 
Berne et al . |2008 1 1 Joblin et al.|2010] >. The region is part of a small 
molecular cloud at 430 pc, illuminated by the Be star HD 200775 
[RA(2000) = 21h01m36.9s ; Dec(2000) = +68° 09 47.8]. It has 
been shaped by the star formation process leading to the for- 
mation of a cavity surrounded by denser filaments and clumps. 
Figure 1 shows the NW part of this region. 



Parameters 




Values 


Units 


Star surface temp. 


T e ff 


15000 


K 


UV radiation field 


Go 


2600 


Habing' 


Fixed pressure 


P 


7x 10 6 


KcirT 3 


A 1 T7 

A V /-CB-V 


R\ 


5.56 






Xo 


4.60 




Extinction curve with 


y 


1.36 


/zrrT 1 


Fitzpatrick & Massa 


C\ 


0.80 




parameters 


C2 


0.32 


yum 


for HD 200775 


C3 


3.09 


/zrrT 2 




c 4 


0.37 




CR ionization rate 


c 


5 x 10" 17 


s- 1 


Dust minimum radius 


^min 


3 x 10~ 7 


cm 


Dust maximum radius 


^max 


3 x 10~ 5 


cm 


MRN distribution index 


a 


3.5 





Notes, t The fluxes are normalized by the flux o f the interstella r radi- 
ation field between 912 and 2400 A measured by Habing ( 1968 1 in the 
~ 1 kpc neighbourhood of the Sun. 




36.00s 34.00s 32.00s 

RA (J2000) 



30.00s 21h01m28.00s 



Fig. 1. Map of the NW PDR of NGC 7023 observed by Spitzer- 
IRAC at 8 yum. The white star shows the location of the exciting star 
HD 200775. The white line shows the Star-NW cut that is studied in 
this paper, and the white crosses indicate the positions at 10, 42 and 
80" from the star. 



|Pilleri et al.| ( (20T2] > used their spectral fitting tool PAHTAT 
to analyse, in this region, the AIB emission that was recorded 
by the Infrared Spectrograph (IRS) on-board Spitzer. They de- 
rived the respective contributions of PAH cations (PAH + ), neu- 
tral PAHs (PAH ) and evaporating very small grains (eVSG). 
Figure [2] reports these values along the Star-NW cut shown in 
Fig. [T] In the following, we focus on two regions of the neb- 
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Table 2. Local physical conditions in the cavity of NGC 7023. 



d 


T 


n(H) 


n(H 2 ) 


n(e-) 


["] 


[K] 


[ cm- 3 ] 


[ cm- 3 ] 


[ cm" 3 ] 


10 


103 


150 


6.9 x 10" 8 


3.7 x lO" 2 


20 


129 


150 


6.8 x 10- 7 


3.9 x lO" 2 


30 


150 


150 


3.3 x 10" 6 


4.0 x lO" 2 



ula along the Star-NW cut: the cavity, where the emission of 
fullerene C(,o is detected and the AIB signal disappears (Sell 



|gren e t al. 2010} |Berne &~Ti elens 2012) and the well-studied 
NW PDR, where the AIB emission is intense. 

In order to provide a realistic description of the spatial evo- 
lution of physical conditions (density and gas temperature) and 
chemical abundances in the PDR, we used the Meudon PDR 
code (version 1.4.3, rev. 608, |Le Bourlot et al.||20L2] ). The 
input parameters are summarized in Tab. [1J The illuminating 
star was modelled using a synthetic stellar atmosphere spectrum 
from Castelli & Kurucz| ( |2004[ ) with an effective temperature of 
15000 K (solar elemental abundances and log(g)=4.0) and an in- 
tegrated UV intensity of 2600 in units of Habing ( |Habing| 1968 ) 
at 42" from the star ( Joblin et al.|2010 1. The UV spectrum mea- 
sured with IUE and the photometric U, B, V ( |Mermilliod et al.| 
|1997| i and J ( |Skrutskie et al.|[2006| l bands were combined with 
this stellar atmosphere to determine the Fitzpatrick and Massa 
parameters of the extinction curve along the line of sight toward 
HD 200775 (Fitzpatrick & Massa 2005). These parameters are 
used for radiative transfer calculations in the PDR code. For 
the NW PDR, we computed an isobaric model with a thermal 
pressure of 7 x 10 6 Kcm 3 , which provides a density profile in 
reasonable agreement with the profile derived by |Pilleri et al.| 
(2012) from the Spitzer mid-IR observations, and also with the 
C + emission at 158 fim measured with HlFl/Herschel at this lo- 
cation ( |Jobhn et al.|2 010). For the cavity, we computed isochore 
models for three positions at 10, 20 and 30" from the star. We 
followed Berne & Tielens (2012) and considered a density of hy- 
drogen nuclei «h = 150 cm" 3 and a local UV radiation field that 
is derived by applying a geometrical dilution factor. 

The resulting physical conditions and chemical abundances 
are shown in Fig. [2] as a function of the position along the Star- 
NW cut for the NW PDR. As the extinction is negligible in the 
cavity, the variations of the radiation field are dominated by the 
geometrical dilution which cannot be taken into account using 
the plane-parallel geometry of the Meudon PDR code. There- 
fore, for each position in the cavity, we ran a different model 
and checked that the extinction was negligible along distances 
on the order of 10". In the following, we use the local condi- 
tions obtained at A v =10- 4 in each cavity model, as summarized 
in Tab. [2] and compute the evolution of PAHs in this region. In 
the next section, we present the physical and chemical processes 
that drive the evolution of PAH charge and hydrogenation states. 



3. Physical and chemical processes 

In the following, we focus on the evolution of the charge and 
hydrogenation states of a few PAHs, namely coronene (Nc=24), 
circumcoronene (Nc=54), circumovalene (Nc-66) and circum- 
circumcoronene (Nc-96). These molecules have compact ge- 
ometries and are therefore expected to be the most stable for a 
given size. Furthermore, coronene is the largest PAH species 



whose chemical properties have been studied experimentally 
( |Bierbaum et al.|2011| ). 

We review here our current knowledge on the main processes 
that drive the evolution of these PAHs, and justify our choices 
when facing the lack of available molecular data. We do not 
consider the gain or loss of C-atoms in this section as (i) reac- 
tivity with gas-phase species such as C + and O is expected to be 
much slower than the processes related to charge and hydrogena- 
tion states ( |Le Page et al.|200T) , and (ii) large and compact PAHs 
were shown to lose primarily all their hydrogen atoms before en- 
during C-loss when submitted to photodissociation (Ekern et al. 
[1998] ). 

In the following, we define the hydrogenation state as the 
number A^h of H-atoms bound to the carbon skeleton of the 
molecule. We disentangle between (i) normally hydrogenated 
PAHs, containing the maximum number A 7 ^ of H-atoms bound 
to peripheral C-atoms without breaking aromaticity, (ii) dehy- 
drogenated PAHs that have lost one or several H-atoms com- 
pared to the normally hydrogenated structure (A^h < A 7 ^), and 
(iii) superhydrogenated PAHs with more H-atoms than in the 
normally hydrogenated structure (A^h > A^), leading to some 
aliphatic bonds. In practice, we will consider only species with 
Nh = A^ + 1 since there is no quantitative data on more hydro- 
genated species. 



3.1. Absorption of UV-visible photons 

The UV-visible absorption cross-sections of PAHs containing 
a few tens of C-atoms have been investigated experimentally 
(|Verstraete et al.|[T9"9"0"} |Joblin et al~||1992| |Ruiterkamp et al.| 
|2002[ ). Larger species are difficult to handle experimentally and 
their cross-sections are more easily obtained with theoretical 
methods, like density functional theory (Mall ocTet al.| |2004). 
One general result of these studies is that UV-visible absorption 
cross-sections scale with the number of carbon atoms in PAHs. 
In this work, we use the cross-sections from the theoreti- 



cal spectral database of polycyclic aromatic hydrocarbons (Mai 
|loci et aL]|2007| ), which provides us with a consistent set of 
data. Since there are no calculations for circumcircumcoronene 
(C96H24), we scaled, according to the number of C-atoms, the 
absorption cross-section of circumcoronene (C54H18), the largest 
studied PAH with similar symmetry. 

3.2. Photoionization 

The photoionization cross-sections of pyrene (CigHio) and 
coronene (C24H12) were measured by |Verstraete et aL ( 1990 », 
who found that they are well parametrised by the number of 
C-atoms and the ionization potential (IP). We make use of the 
formula proposed by Le Page et al. ( |2001| ): 



C~inn — 



0"UV x F ion = cr uv X a x exp{-£>[c(£exc - d)] 4 } 



(1) 



where the UV-visible absorption cross-section cr uv is the one 
discussed in the previous section and accounts for the evo- 
lution with the number of C-atoms. The ionization yield 
Fj on accounts for the effect of the IP. The constants are taken 
from |Le Page et al.| ( |200T| >: a = 0.8, b = 0.00128, c = 
(14.89-IP corone ne)/(14.89-IPpAH) where the ionization potential 
of coronene (IP CO ronene) and of the studied PAH (IPpah) are ex- 
pressed in eV, and d= 14.89 eV. For coronene, circumcoronene 
and circumovalene, we used the IP computed by |Malloci et al.| 
(2007 ) (respectively 7.02, 6.14 and 5.71 eV). For circumcircum- 
coronene, we assumed an IP of 5.68 eV, consistent with the 
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Fig. 2. Left: Analysis of the AIB emission in NGC 7023 along the Star-NW cut, using the PAHTAT tool ( jPilleri et al.|2012| >. The integrated 
intensity over the 5.5 - 14yum spectral range after correction from extinction (thick black), as well as the contributions from eVSGs (red), neutral 
PAHs (green), and PAH cations (blue) are shown. The vertical line indicates the position of the PDR front, at 42" from the star. Right: Physical 
and chemical structure of the NW PDR as computed with the Meudon PDR code, assuming an intensity of the UV radiation field Go=2600 in 
Habing units at A v =0 (42" from the star), and a constant pressure of 7 x 10 6 Kcirr 3 . 



parametrisation proposed by Ruiterkamp et al. (2005) in then- 
Tab. B.l for molecules with a Dg/, geometry. 



3.3. Photodissociation 



As discussed by Allain et al. ( 1996b i, after the absorption of a 
UV-visible photon, a PAH can undergo several transformations 
involving a wealth of processes: it can lose an electron or a frag- 
ment (H, H2, C2H2, ...), or relax energy by radiating IR pho- 
tons or by luminescence. | Allain et al.| ( |1996b[ ) showed that the 
branching ratios between the different pathways vary with the 
internal energy and the size of PAHs. For compact PAHs as 
studied here, the main processes apart from photoionization, are 
H-loss and IR photon emission. In the following, we therefore 
describe the photophysics of PAHs with a few simplified steps: 
(i) absorption of a UV-visible photon of energy hv either leads to 
ionization with the branching ratio Y loa (hv), or increases, by the 
energy hv, the internal energy of the PAH, which is then found 
in a vibrationally excited state of the electronic ground state; (ii) 
emission of IR photons at the rate k\^(E), with E the internal en- 
ergy of the PAH; (iii) fragmentation by H-atom lost at the rate 
kdiss(E). In this section, we present how we determined these 
two rates. 



3.3.1. IR relaxation 

Several approaches have been proposed to compute the IR pho- 
ton emission rat e of PAHs, from the simp le assumption that 



k iR w 10 2 s -1 ( |Herbst & McEwan||1990| l, up to a detailed 



Monte-Carlo kinetics model using a microcanonical formalism 
by |Joblin et al.| ( |2002| l. The latter approach requires the full 
knowledge of the vibrational mode frequencies and associated 
Einstein coefficients. The values of kiz(E) are obtained by av- 
eraging over typically 10000 trajectories. For coronene, circum- 
coronene and circumovalene, we used this method and the vibra- 
tional properties from the theoretical spectral database of PAHs 
(Mall oci et al.|2 007). As can be seen on Fig. [3] the dependence 
of &ir on PAH size is not strong. For circumcircumcoronene, we 
used the same values of Ic\r as for circumcoronene. 



In addition to the determination of kjn(E), we extended the 
method of |Joblin et al.| (2002) to derive the mean energy e^E) 
of the IR photons emitted as a function of the internal energy of 
PAHs, so that one can compute the IR cooling rate as a function 
of the internal energy of the molecule as: 



dE , 
dt 



(E) = k lR (E) x e m (E) 



(2) 



Since the values and shape of k]R(E) were shown to very 
weakly affect the results of PAH evolution modelling (Le Page 
et al. 2001), both values ki^(E) and e^E) were computed only 
for cations in normal hydrogenation state, and used for all charge 
and hydrogenation states. 

3.3.2. Dissociation 

Photodissociation of PAHs has been experimentally and theo- 
retically studied by several teams. Ekern et al. (1998) quali- 
tatively studied PAH photostability and showed that larger and 
more compact species are more resistant to UV-visible photons, 
in agreement with previous results (Jochir ris et al.|[T994] ). The 
same team determined the branching ratios of the dissociation 
paths of fluorene cation (Ci3Hiq + ) showing that the successive 
loss of H-atoms widely dominates the other paths like H2-I0SS 
or C2H2-I 0SS (fD ibben et al.|2001| |Szczepan ski et al.|2001] >. The 
works of|Jochims et al.| ( |1994] i, [Ekern et al.| ( |1998| l and jBaii 1 ] 



isaukas et al.| (|2004) strongly suggest that this behaviour is even 
more pronounced for larger PAHs. The photodissociation rates 
of naphthalene (CioHg + ), anthracene (C^H^), phenanthrene 
(Ci4Hio + ) and pyrene (Ci6Hio + ) have been quantitatively de- 
termined by Lifshitz and co-workers ( |Ho et alT]|1995[ |Ling &| 
|Lifshitz|1998[ [Ling et al.|1995) > who combined experimental re- 
sults with the Rice-Ramsperger-Kassel-Marcus (RRKM) statis- 
tical theory. They showed, in addition, that once a first H-atom 
is lost, the loss of a second H-atom is faste r (|Ling & Lifshitz 
|1998| l. Based on these results, Le Page et al. (2001 1 used a sim- 
plified version of RRKM in order to compute the dissociation 
rates of large PAHs, bearing up to 200 C-atoms. 

The coronene cation (C24Hi2 + ) has been more specifically 
studied by |Joblin et al.| ( |2013b, in prep.[ l with the Fourier 
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Fig. 3. Rates for IR emission and dissociation as a function of the internal energy for coronene C24H12 (left) and circumovalene C66H20 (right). 
For the typical energies of UV photons in PDRs, our results predict dissociation rates that are much lower than in |Berne & Tie lens (2012;), and 
much higher than in Le Page et al. ( 200TJ. 



transform-ion cyclotron resonance trap (FT-ICR) experiment 
PIRENEA. They confirmed that full dehydrogenation occurs be- 
fore any carbon loss. In addition, the partially dehydrogenated 
cations bearing an odd number of H-atoms are found to pho- 
todissociate much faster than those bearing an even number of 
H-atoms, consistently with the results of Ling & Lifshitz ( 1998) 
on smaller species. 



The results of Jochims et al. ( 1994[ > and Ling & Lifshitz 



1998) are well explained by statistical theories. This justifies 



the use of statistical approaches in previous PAH models (Leger 



et al. 



1989; |Allain et al.|1996b||Le Page et al.|2001||Visser et al 



2007 1. We chose to follow the approach described in|Leger et al. 



1989 ) and Boissel et al.| ( |1997] ), based on the Laplace transform 



of the Arrhenius law: 



kdiss(E) = A d 



pCg - £0) 

S P(E) 



(3) 



where E is the internal energy of the PAH, Ad; ss has the dimen- 
sion of a frequency, Eq is the dissociation energy and p(E) is 
the vibrational density of states of the PAH. This very simple 
approach was shown to provide satisfying results when com- 
pared to more sophisticated methods, like RRKM theory ( |Barker| 
[1983) . 

The two parameters Adi ss and Eq were discussed for the dif- 



ferent hydrogenation states of coronene by Jobli n et aL] ( 2013b,| 
in prep. ). The authors show that assuming a dissociation en- 



ergy £o=4.8 eV for the C-H bond of all PAHs bearing an even 
number of H-atoms, they manage to account for the experi- 
mental results for all these species by using a single value of 
Adiss = 6.8 x 10 17 s -1 that was derived from the photoion appear- 
ance potential measured by |Jochims et al.| f!994) on €24^2. In 
addition, the dissociation of PAHs bearing an odd number of H- 
atoms is well described with the same constant A d i ss by assum- 
ing Eo = 3.2 eV. Since there is no data available for the other 
species, we use, in the following, the same values of Adi ss and 
Eo for the four PAHs considered in this work, in their cationic 
and neutral states. Therefore, the evolution of fcdiss with molec- 
ular size is driven by the dependence of the vibrational densities 
of states on size, which basically reflects the evolution of the 
number of degrees of freedom. We computed the exact quan- 
tum harmonic vibrational densities of states with the method of 
Beyer & Swinehart ( jStein & Rabinovitch|1973] ). 



Data concerning the dissociation of superhydrogenated 
PAHs are scarce. They are thought to be more easily disso- 
ciated than molecules in a lower hydrogenation state (Le Page 
|et al.|200"T] and references therein). Therefore, in the following, 
we use the dissociation rate of the corresponding PAH missing 
a single H-atom as a lower limit for their dissociation rate (e.g. 
£diss(C24Hi3) = &diss(C24Hn) ), leading to an upper limit on the 
abundances of superhydrogenated species in our results. 

Compared to other studies, our work provides significantly 
different values. The new experimental constraints on coronene 
cations lead to much more fragile species than predicted by 
|Le Page et al.| ( |2001| ) (see Fig. [3). On the opposite, the more 
empirical method that is presented in |Tielens (2005 1 and used in 
[Berne & Tielens| ( f2012| ), leads to very different values and shapes 
of £diss(£)> and a lower stability compared to our values. We ex- 
pect our dissociation rates to provide better results because they 
are based on experimental data on larger PAHs than considered 
in previous studies. 



3.4. Reactivity with hydrogen 

Only few experimental studies on the reactivity of gas-phase 
PAHs with hydrogen have been published, and all of them con- 
cern small cations bearing up to 24 C-atoms. Therefore, we eval- 
uate first the reaction rates for coronene cations in all their hy- 
drogenation states from the available data, and then extrapolate 
these values for the larger PAHs considered in our work. The 
values of reaction rates for C24Hn + and their origin are summa- 
rized in Tab. [3] 

Reactions between PAHs and atomic or molecular hydro- 
gen generally lead to the addition of one or two H-atoms to pe- 
ripheral C-atoms of the PAH (see discussion by Le Page et al. 



2001). The reactivity of PAH cations with atomic hydrogen 



is found to be rather fast and to decrease when PAH size in- 



creases: 



1992), 1.9 x 10 



- 10 cm 3 s- 1 



for benzene (CgH6 
for naphthalene (C10H 
1 for pyrene (Ci6H 10 + ,|Le Page et al. 



Petrie et al. 
iSnow et al. 



for coronene (C24H17T Betts 



M998]) 1.4 x 10" 1U cm 
1999a|), an d 1.4 x 10" 
|et al.|2006) . 

Following Le Pag e et aT] ( |2001[ ) we use, for normally hy- 
drogenated PAH cation C x H y + , a reaction rate that scales with 
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Table 3. Summary of the measured and inferred reaction rates of atomic and molecular hydrogen with coronene cations in their various hydro- 
genation states. 



Reaction 


k [cm 3 s 


] 


Ref. 


Comments 




Evolution with Nc 


C24H12 + H 


1.4 x 10~ 


10 


(a) 


FA-SIFT 




It 1 X 


C24Hi2-2n + + H 


1.4x10" 


10 


(b,c,d) 


extrapolation from CioH6 + 


and C 16 H 8 + 


y/x 


C24Hi2-2n+l + + H 


~ 5 x 10- 


1 1 


(b,c) 


extrapolation from CioH7 + 




constant 


C24Hi2 +n + + H 


~ 10- 12 




(e) 


extrapolation from C(,Hj + , 


Ci H 9 + and Ci 6 H n + 


constant 


C24Hi2 + + H2 


< 5 x 10- 


u 


(a) 


FA-SIFT 




constant 


C24Hi2-2n + + H2 


< 5 x 10- 


13 


(b) 


extrapolation from CiqH6 + 


and C 16 H 8 + 


constant 


C24Hi2-2n+l + + H2 
C24Hi2 +n + + H 2 















Notes. The type of experimental set-up used to perform the measurements is given in the "Comments" column. When "extrapolation" is mentioned, 
the rate value was extrapolated from measurements on smaller species. FA-SIFT: flowing afterglow-selected ion flow tube. The last column 
mention the kind of variation of the rates with the size of PAHs that we considered in our model. When y/x is mentioned, the rate varies 
proportionately to the ratio of peripheral (y) to total (x) C-atoms in the species. 

References. (a)|Betts et ah]([2006li; (b)|Le Page etal^pOOT); (c)|LePage etaL](|1999a[>; (d)|Le Page etaLl<|1999b|>; (e)|Snow et afT]fT998>. 



^coronene ■ 



the reaction rate of coronene, according to k VJ - (y/x)k, 
This expression is in reasonable agreement with available exper- 
imental results considering a typical experimental uncertainty 
of 50% (|Betts et al.|[2006]>. We also follow |Le Page et"aL 



( |2001| l to estimate reaction rates for dehydrogenated and super- 
hydrogenated PAHs, assuming that (i) dehydrogenated coronene 
cations missing an even number of H-atoms react with the same 
rate than the normally hydrogenated coronene cation, as inferred 
from the behaviour of naphthalene and pyrene dehydrogenated 
cations; (ii) dehydrogenated coronene cations missing an odd 
number of H-atoms react with a lower rate of ~ 5 x 10"' 1 cm 3 s -1 
similarly to CioH7 + ; (iii) super-hydrogenated coronene has a 
lower rate of ~ 10~ 12 cm 3 s- 1 , as suggested by measurements 
on smaller PAHs (Sno w et al.||1998"] l. The reaction rates for 
the other PAH species are derived by using a similar scaling for 
even-de hydrogenated PAHs as for normally hydrogenated PAHs. 
No evolution of reaction rates with size is considered for odd- 
dehydrogenated PAHs and superhydrogenated PAHs. 



The reaction rates of PAH cations with molecular hydrogen 
were determined to be below the experimental detection thresh 
old, which led to an upper limit of 
reaction rate of C24H12 



5 x 10 13 cm 3 s" 1 for the 



+ H 2 (|Betts et al. 2006 1. We there 



fore calculate upper and lower limits to the abundances of PAH 
hydrogenation states by assuming either that PAHs do not react 
with H2, or that PAH cations react with H2 with a reaction rate 
of 5 x 10 13 cm 3 s -1 leading to the addition of two H-atoms to 
peripheral sites. 



3.5. Recombination with electrons 

Recombination of PAH cations with electrons leads to an in- 
crease of the PAH internal energy by the value of the ionization 
potential. The statistical approach presented in Sect. 3.3.2 can 



then be applied to determine whether this extra energy is relaxed 
radiatively or triggers molecular dissociation. According to the 
rates shown in Fig. [3] in the case of coronene, IR photon emis- 
sion is at least 100 times faster than dissociation for an internal 
energy of 7 eV, and the recombination of coronene cations with 
electrons is therefore expected to be non-dissociative. Larger 
species have smaller values of their ionization potentials and are 
more stable than coronene. We can therefore neglect the disso- 
ciative channel for all the species studied here. 

The recombination rates of PAH cations with electrons were 
measured for a numb er of small PAHs ([Abouelaziz et al.| 1993 



|Hassouna et al. 2003 ; R ebrion-Rowe et aL[2 003 ; Novotny et al. 
| 2005| |Biennier et al. 2006 1. All the experimental values are 



about one order of magnitude lower than the Spitzer law, which 
is based on the classical model of a thin conducting disk ( Spitzer 
|1978|fVerstraete et al.|19 90). However this difference decreases 
when the size of PAHs increases. To our knowledge, there is 
no published model available that accurately accounts for this 
trend. Therefore, for the four PAHs of this study, we assume that 
the recombination rates are comprised between the experimental 
rate of pyrene cations and the Spitzer law, and use the average of 
these two values in our model. This leads to rates roughly one 
order of magnitude higher than those proposed by Le Page et alT] 
( 12001) . 



The reactivity of neutral PAHs with atomic hydrogen is un- 
known. Mebel et al. ( 1997| l theoretically showed the absence of 
a potential barrier for the reaction CgH6 + H, suggesting that re- 
action rates may be high. However, Le Pa ge et al7] ( |200 1 ) > found 
that the results were not affected when including a reaction rate 
of neutral PAHs with atomic hydrogen of 10 cm 3 s _1 . In our 
standard model, we neglect the reactivity of neutral PAHs, and 
evaluate afterwards the uncertainty resulting from this choice, 
arbitrarily assuming reaction rates of 3 x 10" 11 cm 3 s -1 at 300 K 
for radical neutral PAHs (odd number of H-atoms), and 3 x 10" 13 
cm 3 s 1 at 300 K for closed-shell neutral PAHs (even number of 
H-atoms). 



3.6. Charge exchange 

Charge exchange reactions between neutral PAHs and ions such 
as C + , Fe + or Si + may impact the ionization bal ance of PAHs. 
Reactions with C + proceed at the Langevin rate ( Canosa et al. 
1995), and can marginally affect the PAH charge state in low 
UV irradiation environments like the diffuse ISM (Wolfire et al. 



2008). Photoionization rapidly dominates the ionization balance 
of PAHs when the UV radiation field intensity increases. Other 
possible partners for charge exchange are much less abundant 
than C + and are therefore completely negligible. In the follow- 
ing, we do not consider charge exchange reactions. We will 
check the validity of this assumption in Sect. |5.3| 
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Table 4. Species and main parameters of the processes considered in our standard model. 



N C K q N H 


IP a (300 K) 
[eV] [cmV] 


Eq ^diss 

[eV] [s- 1 ] 


fc +H (300K) k +il2 
[cm 3 s~ I ] [cm 3 s~'] 




In 

2n + 1 

& 24 12 13 
" 2n 

+ 2n + 1 

13 


7.02 
7.02 
7.02 

1.0 xl0~ 5 
1.0 xl0~ 5 
1.0 xlO" 5 


4.8 
3.2 

32 6.8 x 10 17 
4.8 

3.2 

3.2 


0/3 xlO- 13 
0/3x10-" 
0/3 xlO- 13 
1.4 xlO" 10 0/5 xl0~ 13 
5x10"" 0/5x10-° 




In 

jju 2« + 1 

& 54 18 19 
?gC 2n 

+ 2n + 1 

19 


6.14 
6.14 
6.14 

1.4 xl0~ 5 
1.4 xlO" 5 
1.4 xl0~ 5 


4.8 
3.2 

32 6.8 x 10 17 
4.8 

3.2 

3.2 


0/3 xlO- 13 
0/3x10-" 
0/3 xlO- 13 
9.3x10-" 0/5X10" 13 
5x10-" 0/5 xlO- 13 




2« 

» 2« + 1 

VtV 66 20 21 
OX 2n 

+ In + 1 

21 


5.71 
5.71 
5.71 

1.56 x 10~ 5 
1.56 xlO" 5 
1.56 xlO" 5 


4.8 
3.2 

32 6.8 x 10 17 
4.8 

3.2 

3.2 


0/3 xlO- 13 
0/3x10-" 
0/3 xlO" 13 
8.5x10-" 0/5 xlO- 13 
5x10-" 0/5 xlO" 13 


Wj 


2« 

ftV 2« + 1 

ttft 96 24 25 
IIC 2n 

W + 2« + 1 

25 


5.68 
5.68 
5.68 

1.85 x 10~ 5 
1.85 xlO" 5 
1.85 x 10~ 5 


4.8 
3.2 

32 6.8 x 10 17 
4.8 

3.2 

3.2 


0/3 xlO" 13 
0/3x10-" 
0/3 xlO" 13 
7.0x10-" 0/5 xlO" 13 
5x10-" 0/5 xlO" 13 



Notes. The columns give: (1) the structure of the species studied in this paper in their neutral normally hydrogenated state; (2) the number Nq of 
C-atoms in the species; (3) the number of H-atoms in the normally hydrogenated state; (4) the charge q of the species; (5) the number N a of 
H-atoms in the species; 2n is for normally hydrogenated PAHs and PAHs missing an even number of H-atoms, 2n + 1 is for PAHs missing an even 
number of H-atoms; other values are superhydrogenated states; (6) the ionization potential IP of the species; (7) the rate of recombination a with 
electrons at 300 K; (8) the binding energy E and (9) the pre-exponential factor A diss used in Eq.|3] (10) the reaction rate for H-atom addition at 
300 K; (11) the rate for reaction with H 2 . For (10) and (11), when two values are given, they correspond to lower/upper limits considered in our 
model. 



3. 7. Sum-up of our standard model 

The species and processes included in our model are summa- 
rized in Tab. |4] In short, we consider here four sizes of PAHs 
in various hydrogenation states and two charge states, namely 
coronene (C24Hn°^ + ), circumcoronene (C54Hn°^ + ), circum- 
ovalene (C66Hn°^ + ) and circumcircumcoronene (C96Hn°^ + ). We 
limit the possible charge states to neutrals and monocations, ex- 
cluding dications and anions for several converging reasons, in- 
cluding (i) the lack of data for the recombination of dications 
with electrons and for the association of neutral PAHs with elec- 
trons, and (ii) the expected weak abundances of these other 
charge states in the environments considered here. The latter 
point can be checked afterwards. We consider all the hydrogena- 
tion states from fully dehydrogenated species, i.e. pure carbon 
clusters, to superhydrogenated PAHs bearing one extra H-atom. 
This latter choice is motivated by the lack of quantitative data 
for both the dissociation rates of superhydrogenated species and 
the rates of reaction of these species with hydrogen. The disso- 
ciation rate that is used for superhydrogenated PAHs should be 
considered as a lower limit. 

The charge evolution is driven by the balance between pho- 
toionization and recombination with electrons, both depending 



exclusively on the size of PAHs, in our model. Similarly, the 
hydrogenation state is driven by the balance between photodis- 
sociation and reactivity with atomic and molecular hydrogen. 
Photodissociation strongly depends on the size of the species 
through their vibrational density of states, and also strongly on 
their hydrogenation state through the binding energy Eq, PAHs 
with an odd number of H being less stable than PAHs with an 
even number of H. Dependence on the charge state is not taken 
into account in our model for this process. On the contrary, re- 
activity with atomic and molecular hydrogen mainly depends on 
the charge and, to a lower extent, on the hydrogenation state of 
PAHs, and more marginally on their size. Considering the lack 
of data for the reaction rates of PAH cations with H2 and of neu- 
tral PAHs with H, we consider lower and upper values for these 
processes. 

Describing PAH evolution requires to model the processes 
discussed above, as well as how they vary with astrophysical 
conditions. The UV-visible radiation field spectrum drives both 
photoionization and photodissociation. The abundances of H, 
H2 and free electrons are involved through their reactions with 
PAHs, in which the gas kinetic temperature T plays a role. Con- 
sidering only the influence of temperature on the collision rates, 
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the rates of recombination with electrons vary like T~ l l 2 (see, 
e.g., Bakes & Tielens 1994), the rates of reaction between two 
neutral species like T /z , and no variation with T is considered 
for reactions between charged and neutral species (Le Page et al. 
|200Tj |Herbst|200i1 and references therein). 



4. Numerical framework 

4.1. Description of the species 

The species considered in this work are presented in Sect. |3.7| 
and Tab. [4] The evolution of these species involves some pro- 
cesses that critically depend on the internal energy of the re- 
actants, like photodissociation processes. An explicit descrip- 
tion of this internal energy was proposed by Visser et al. ( 2007 ) 
in the context of protoplanetary disks, where the strong UV- 
visible irradiation by the close illuminating star is expected to 
favour multiphoton events, i.e. absorption of another photon 
before complete cooling of the PAHs. Previous studies on the 
evolution of PAHs in regions of low or moderate excitation do 
not describe explicitly the evolution of the internal energy of 
PAHs (|Allain et al.|[l996al |Le Page et aLl[200T] |Berne & Tie-| 
|lens|2012| >. Nevertheless, a quick examination of the rates for IR 
photon emission (&ir) and dissociation (&diss> see Fig- [3]) reveals 
that for large enough species, the threshold energy E& for which 
^iR(^th) - ^diss(^th) falls we U above the Lyman cut-off at 13.6 
eV. As a consequence, even in low radiation fields, dissociation 
is driven by multiphoton events, since this is the only way to 
raise the internal energy of large PAHs close to or above E^. 

In order to describe the evolution of internal energy, each 
species X was divided into an ensemble of subspecies X, (i=0, 
. . . , iVbin) whose number density is the number per unit of vol- 
ume of species X, bearing an internal energy between (i - 1) X AE 
and i x AE for (i=l, . . . , Nbm) or no internal energy for i - 
0. Thus, instead of studying the evolution of an ensemble of 
species, we study the evolution of the histograms of their in- 
ternal energy. Obviously, the accuracy of the method depends 
on the width of the bins, AE, that we found to be optimum for 
AE = 0.25 eV, which is a compromise between accuracy and 
computing time (see appendix [A4]i. 



the following scheme: 

„ , [l-JW0]*«b 5 (0 

Xi + hv m — > Xy 



4.2. Evolutionary scheme 

We aim at computing the time evolution of the species defined in 
the previous section in a specific environment. The feedback of 
PAH evolution on their astrophysical environment is out of the 
scope of this work. Therefore, the densities of H, H2 and free 
electron are considered as fixed external parameters, as are the 
kinetic gas temperature and the UV-visible flux. We compute 
the time evolution of each species X and its internal energy his- 
tograms by considering the chemical balance of each subspecies 
X t : 



d[Xj] 
dt 



= Pi - h 



(4) 



where [X,] is the number density of the species X,, P, and L, 
are the production and loss rates, respectively. These two latter 
quantities are computed from (i) the conditions of the astrophys- 
ical environment, and (ii) the rates of the processes presented 
in Sect. [3] Our work differs from previous studies by the way 
we model photodissociation processes. Instead of computing an 
effective photodissociation rate for each species X, we compute 
the production and loss rates for each subspecies Xi according to 



*r(0 



X v , + hv.. 



(5) 
(6) 
(7) 



where i' > i, i" < i, and X Q H has no internal energy. The nu- 
merical rates [1 - Tj on (0]£ a bsOX &ir(0 and fcdissO) are computed 
from the physical quantities Y\ on (E), cr a b s (£'), the UV-visible 
flux, k 1R (E) and k^ ss (E), as detailed in appendix[A] For the other 
processes, we do not take into account the effects of the inter- 
nal energy of PAHs, which makes the numerical implementation 
straightforward. 

5. Results 

The computed charge and hydrogenation states of coronene, cir- 
cumcoronene, circumovalene and circumcircumcoronene in the 
NGC 7023 NW PDR are shown in Figs. H] and [5]for different as- 
sumptions concerning the reactivity of PAHs with hydrogen. For 
all our calculations, we assumed neutral normally hydrogenated 
PAHs as initial conditions. In this section, we present our results 
focusing successively on the charge state and the hydrogenation 
state of PAHs once steady state is reached. 

5.1. Evolution of PAH charge in NGC 7023-NW 

Using our model, we computed the global cation fraction for the 
four PAHs considered in this study by summing the cation frac- 
tion in all hydrogenation states (e.g. for coronene £jv h C24Hn h + 

/ 2a- h , 9 Qz4H Nh « with n h = P. 12 D- m Fi S- the calculated 
values are compared with the cation fraction extracted from the 
analysis of the AIB spectrum along the Star-NW cut (Sect. |2]i. 

The computed cation fractions are in good agreement with 
the observed one in the cavity, as we predict fractions above 
90% of cations at 20" from the star. In the region of the PDR 
where ionization of PAHs is found to be significant (Ay<; 2), 
photoionization dominates the evolution of the charge state of 
PAHs since the UV radiation field intensity varies by a factor 
~ 10, whereas the PDR model predicts a rather constant abun- 
dance of free electrons and a gas temperature that changes by a 
factor of 2 (see Fig. [2}, which results in a variation of the rate 
of recombination of PAH cations with free electrons by a factor 
of less than 2. The agreement between computed and observed 
cation fractions is also good within the AIB peak emission region 
(between 42" and 52" from the star along the Star-NW cut), both 
in terms of absolute values and of profile evolution. Finally, the 
cation fraction is found to increase with the PAH size, Nc = 96 
being twice more ionized than A^c = 24. A more quantitative 
comparison between modelling and observations would require 
to weight each calculated charge ratios by the column density of 
the corresponding PAH along the line of sight. 

5.2. Hydrogenation state 

As shown in Fig. |4] the hydrogenation state of PAHs evolves 
significantly with both PAH size and position within the PDR. 
A general result is that partially dehydrogenated species never 
dominate in terms of abundances. Even collectively, their maxi- 
mum weight is 40%, and for a very narrow range of spatial posi- 
tions. Therefore, the resulting hydrogenation states can be clas- 
sified into four categories: pure carbon clusters (A^h = 0), par- 
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Fig. 4. Charge and hydrogenation states (number of H-atoms in the PAH) of coronene (C24H12), circumcoronene (C54H18), circumovalene 
(C 66 H 2 o) and circumcircumcoronene (C96H24) as a function of the position in the NW PDR of NGC 7023, with the astrophysical conditions 
presented in Fig.|2]and with fc+n, = and A: + H(neutral) = 0. Neutral species are in red, cations in blue. 
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Fig. 6. Fraction of ionized species (PAHs and carbon clusters) along 
the Star-NW cut in the NW PDR as derived from the observations (yel- 
low, from[Pilleri et al. 2012} and as computed with our model (other 
colours). The total AIB flux is also plotted in black. The grey surface 
shows the uncertainty range on the ionized fraction for Cs4Hn + species 
(N=0-19), which reflects the range of values we considered for the elec- 
tron recombination rate. The PDR front is assumed to lay at 42" from 
the star. 



tially dehydrogenated PAHs (1 < Nn < A 7 ^ - 1), normally hy- 
drogenated PAHs, and superhydrogenated PAHs (Nn = iV^+1). 
This behaviour results from a general trend that more dehydro- 
genated PAHs are less photostable. Therefore, once a normally 
hydrogenated PAH starts losing H-atoms, further dehydrogena- 
tion is even faster. 

We summarize the evolution of the dominant hydrogenation 
state as a function of the position and size of the species in Fig. [7] 
(results in the cavity are not shown). To illustrate the contribu- 
tion of each population to the AIB emission, Fig. [7] also shows 
Paibi the emission power per unit volume in the mid-IR domain. 
It was calculated using the formula Paib = £"h Go ( |Pilleri et! 
2012 1, where e = 5 x 1(T 32 W H 1 and Go is expressed in units of 
Habing. Note that in this formula, all species, whether PAHs or 
carbon clusters, are assumed to emit the same power per C-atom. 
For a given position in the PDR, larger molecules tend to be more 
hydrogenated than smaller, as expected. For a given size, how- 
ever, one find alternatively hydrogenated and dehydrogenated 
species when increasing the distance from the star, whereas a 
monotonous behaviour, as a function of the intensity of the UV 
radiation field, could be expected. Similarly, superhydrogenated 
PAHs dominate close to the cloud surface, and are less abundant 
in deeper layers. These inversion features are due to the decrease 
in atomic hydrogen density from Av=0.2 to Ay~ 2, and is there- 
fore a consequence of both the environment evolution and the 
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Fig. 5. Same as Fig.|4]but with a reactivity of PAH cations with IT of k + u 2 = 5 x 1CT 13 cm 3 s -1 (upper panel), and with a reactivity of neutral 
PAHs with H of A:™^ = 3 X 10~ u cm 3 s _1 for radicals and fe^ = 3 X 10~ 13 cm 3 s -1 for closed-shell (lower panel). Figures for coronene (C94Hn +/0 ) 
and circumcircumcoronene (C96H N +/0 ) are not shown as they exhibit only marginal variations with respect to Fig. [4] Neutral species are in red, 
cations in blue. 



assumption that reactivity with H2 can be neglected. We discuss 
this assumption further in Sect. 



6.1 



Superhydrogenated PAHs are found to be abundant only for 
the largest species (Nc = 96), despite the fact that we con- 
sidered a lower limit for the dissociation rate of superhydro- 
genated PAHs. Pure carbon clusters are found to be the dom- 
inant form for small- (Nc < 50) and, possibly, medium-sized 
(50 < Nq ^ 70) species, until cloud depths, at which PAHs are 
thought to be part of very small grains (Ay>2, see Fig. [2] and 
Pilleri et al.|2Q12) >. In the cavity, where the density is lower and 
the radiation field much stronger, we predict all PAHs to be com- 
pletely dehydrogenated. 

Interestingly, circumcoronene and circumovalene are found 
to be fully dehydrogenated in the NW PDR at Av=0, and be- 
tween Av=l and Ay =2, respectively, whereas for the same con- 
ditions previous models ( |Le Page et al.||200Tj |Berne & Tielens| 
2012} would predict normally hydrogenated circumovalene to 
dominate. This discrepancy is due to multiphoton events that 
are taken into account in our model. Multiphoton events are by 
far the main channel for PAH dissociation at threshold energy 
(where IR photon emission and dissociation are equally proba- 
ble) above the Lyman limit. Figure [8] shows in the case of cir- 
cumovalene C66H20 that cooling by IR photon emission is 10 
orders of magnitude faster than dissociation at wavelengths of 



maximum UV radiation. C66H20 molecules bearing an internal 
energy close to the threshold energy (~ 20 eV) after the succes- 
sive absorption of several UV photons are found to be the most 
numerous to dissociate, despite their very low abundances (cf. 
Figure [HJ. 

5.3. Timescales 

The results presented in the previous sections are obtained after 
steady state has been reached. We investigate here how steady 
state timescales compare with other evolutionary timescales in a 
given environment. The steady state for PAH charge evolution 
is reached within a few hours at the position A v =l. As a com- 
parison, for these conditions, the typical timescale for C pho- 
toionization is ~1 month fLe Petit et al.||2006| >. This timescale 
can also be compared with the timescale for charge exchange 
(c.e.) between neutral PAHs and C + . In the cond itions of NGC 
7023 N W, using k ce = 3 x 10~ 9 cm 3 s" 1 from |Canosa et al. 
( |l995t , h h = 10 4 cm 3 , C/H = 1.3 x 10 4 from our PDR calcu- 
lations, and assuming that all atomic carbon is ionized, we find 
t = l/(£ c . e .[C + ]) ~ 8 years. This shows that charge exchange is 
negligible in these conditions. 

The steady state for PAH hydrogenation is reached for all 
species on much longer timescales, which span a very broad 
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Fig. 8. Upper panel: Spectrum of the local UV- visible radiation field 
in NGC 7023-NW at 42" (A v =0) from the star, compared to the rates 
of dissociation and IR emission of the circumovalene molecule C66H20- 
Lower panel: The internal energy of C6 6 H 2 o (red histogram) computed 
with our model in steady state, in the radiation field showed in the upper 
panel. The black histogram shows the distribution of internal energies 
of C66H20 molecules prior to dissociation. The vertical dashed lines 
show the 13.6 eV limit, emphasising that dissociation is dominated by 
multiphoton events. 
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Fig. 7. Dominant hydrogenation state as a function of PAH size and 
position in the NGC 7023 NW PDR, without reactivity {upper panel) or 
with a rate of 5 x 10~ 13 cm 3 s~' (middle panel) for the reaction between 
PAH cations and H2. The points represent the parameters for which 
models were computed. Red, yellow and green represent carbon clus- 
ters, normally hydrogenated PAHs, and superhydrogenated PAHs (lim- 
ited to one extra H-atom), respectively. Lines delimiting the different 
domains were interpolated to guide the eye and should be considered 
with caution. The power per unit of volume P Am emitted in the mid-IR 
domain by PAHs and carbon clusters as a function of their position is 
shown in the lower panel, and is also overlaid as a brightness scale to 
the upper and middle panels. The variation observed between Av=0 and 
0.2 reflects the variation of n H ( see Fig. [2]). 



range of values depending on PAH size. Small species reach 
there steady state within less than a year (e.g. ~1 month 
for coronene at Ay-l), while large species exhibit very long 
timescales, as can be seen in figure|9] with ~ 10 2 years at Ay=0 
and ~ 10 4 years at Ay=l in the case of circumovalene, if no re- 
activity is considered with H2. When considering the reactivity 
of PAH cations with H 2 at a rate of k + n, - 5 x 10" 13 cm 3 s _1 , 
the timescale at Ay- 1 decreases down to ~ 10 2 - 10 3 years. 
These values are very large with respect to the typical hydro- 
genation time (l/(k+u x «h) ~ 10 days), and seem surprising 



at first sight. In fact, this is linked to the absence of dominat- 
ing partially hydrogenated species: the global time evolution for 
one PAH species scales with the typical times needed for it to 
evolve through all its partially hydrogenated states, from the nor- 
mally hydrogenated state to the fully dehydrogenated state and 
vice versa. These timescales become long when multiphoton 
events are required to dissociate PAHs, as each dehydrogenation 
step behaves as a limiting step. This illustrates that the complex 
evolution of the hydrogenation state of PAHs requires a detailed 
modelling. 

Considering that |Alecian et al.| ( |2008[ ) estimated the age of 
HD 200775 around lO 3 years, our results suggest that the evolu- 
tion of large PAHs is coupled to the dynamics of the region, and 
that large PAHs may never reach their steady state abundances. 



6. Sensitivity to parameters 

6. 1 . Identifying the key processes 

We have discussed in Sect.[3]the uncertainties on the parameters 
of the various processes. In this section, we aim at evaluating the 
sensitivity of the results to some of them, in order to identify the 
processes that require further investigation. 

The determination of the charge state of PAHs is mainly lim- 
ited by the knowledge of the rate a for their recombination with 
electrons. In the case of circumcoronene, we computed addi- 
tional models with the upper and lower values of a as defined in 
Sect. |3.5| In figure [6| we represented by a grey area the domain 
of values resulting from this uncertainty. Thus, the fraction of 
PAH cations is determined with a typical uncertainty of a fac- 
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Fig. 9. Time evolution of some charge and hydrogenation states of circumovalene in the NGCC 7023 NW PDR, at depths corresponding to A v = 
(a) and Ay= 1 (b), i.e. ~ 42 and ~ 50" from the star, respectively. The typical timescale to reach steady state for charge and hydrogenation states 
are indicated with black arrows. 



tor of two, which is slightly larger than the dispersion of cation 
fraction between PAHs of different sizes. This means that fur- 
ther constraints on the PAH size distribution using PAH charge 
states would require more investigation of the electronic recom- 
bination process. 

The results presented in Fig.|4]do not consider any reactivity 
of PAHs with H2, and of neutral PAHs with H. Figure [5] shows 
the influence of these two processes. The reactivity of neutral 
PAHs impacts the first layers of PDRs, at the interface with the 
cavity where atomic hydrogen dominates, whereas the reactiv- 
ity of PAH cations with H2 plays a significant role deeper in 
the PDR, between A v =0.5 and 3. The smallest in size species 
(Nc = 24) are not affected by these processes as they remain 
fully dehydrogenated. Similarly, the high photostability of the 
largest species (Nc = 96) ensures them to remain normally 
hydrogenated regardless of these additional processes. On the 
contrary, intermediate species can be significantly more hydro- 
genated. The size - position diagrams (Fig. [7]) shows that the 
uncertainty on the reaction rates with hydrogen leads to some 
uncertainty on the position of the transition between normally 
hydrogenated PAHs and carbon clusters, which lays between 
N c ~ 50 and ~ 60 at A v =0 and between N c ~ 50 and ~ 70 
at Ay=l. Therefore, the lack of knowledge on these processes, 
and particularly the reactivity with H2, appears as a major limi- 
tation for predicting the minimum size of PAHs in PDRs. 



6.2. Sensitivity to the astrophysical conditions 



fractions of atomic hydrogen n(H) = 0.5«h and free electrons 
n(0 = 1.4 x 10- 4 n H . 

The results are presented in Fig. [10] They show that the evo- 
lution of the hydrogenation state of C54H18 and C66H20 exhibits 
a complex behaviour with Go and «h- On average, the slope 
of the boundary between the hydrogenated and dehydrogenated 
domains (solid line in Fig. lOl is clearly greater that 1 (~1.5 
and ~2 for circumcoronene and circumovalene, respectively) in 
a log-log diagram. This results from the fact that, on average, 
roughly two photons are needed for C54H1S and C66H20 to reach 
their threshold energies for fragmentation at ~ 17 and ~ 20 eV, 
respectively. When the intensity of the UV radiation field in- 
creases, the proportion of multiple absorption events increases 
and the probability of low energy photons being involved in- 
creases, which leads to slope changes with Go in Fig. [10] 

Coronene is the only species for which we find that the hy- 
drogenation state scales with Go/«h, similarly to the results re- 
ported by Le Pag e et aL] ([2003 ). who did not include multi- 
photon events in their model. This is due to the fact that coronene 
has a small threshold energy for fragmentation (Efh ~ 9.6 < 13.6 
eV) and therefore can easily dissociate after the absorption of a 
single photon. In other words, Go/«h is a correct parameter to 
quantify the hydrogenation state of PAHs only in cases in which 
dissociation is achieved with a single UV photon. However, in 
such cases dehydrogenation is so fast that it leads to the complete 
dehydrogenated product. Therefore, the population of interstel- 
lar PAHs must be dominated by species for which multi-photon 
events dominate the photodissociation process, and whose hy- 
drogenation state is not well parametrised by Go/«h- 



Deriving the general trends of PAH evolution as a function of 
the astrophysical conditions from the case study of NGC 7023- 
NW is not straightforward, because the evolution of the physical 
conditions are correlated along the Star-NW cut. Therefore, in 
addition to the case study of NGC 7023, we also followed the 
method proposed by Le Pa ge et aT] ( |2001| l and computed the 
charge and hydrogenated states of the four studied PAHs on a 
grid of astrophysical conditions. We used the interstellar radia- 
tion field (ISRF) of Mat his et al] ( |1983| l as UV- visible spectrum, 
and normalised it to integrated intensities of 10" (n = 0, 1 , 2, 3, 4) 
in Habing units. We considered a gas temperature of 100 K, and 
several values of «h between 10 and 10 5 cm 4 , with constant 



7. Discussion and perspectives 

7.1. Consequences for the identification of a specific PAH in 
the ISM 

The detection of specific PAHs in the ISM has remained unsuc- 
cessful: |Pilleri e t al. ( 2009 ) reported an upper limit for rotational 
lines of corannulene C20H10 in the Red Rectangle (RR) nebula, 
and |Kokkin et al. (2008) also reported an upper limit on the 
abundance for hexa-peri-hexabenzocoronene C42H18 in the dif- 
fuse ISM from its optical emission. As a comparison, the results 
presented in Fig. [10] show that the larger, more compact, and 
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in Sect. |3.3.T| and computed the dissociation rates for the losses 
of C, C2 and C3 using Eq. [3] and the parameters proposed by 



Leger et al. (T9891 that are based on thermodynamic properties 



0.5 1 1.5 2 2.5 3 
log 10 (G ) [Habing] 



3.5 



Fig. 10. Hydrogenation state of circumcoronene (C54H18, red) and 
circumovalene (C66H20, black) as a function of the intensity of the 
UV radiation field Go in Habing units and the local density «h- The 
lines labelled with 0.1 (respectively 0.5, 0.9) show the points for which 
10% (respectively 50%, 90%) of the species are totally dehydrogenated. 
For coronene, only the blue line labelled with 0.995 where 99.5% of 
coronene is fully dehydrogenated is represented, which is the dominant 
state of coronene for all the conditions spanned here. 



therefore more photostable circumcoronene molecule Cs^^ is 
expected to be at the edge between normally hydrogenated and 
fully dehydrogenated in physical conditions typical of the dif- 
fuse ISM (e.g. Go ~ 1 and n H ~10-100 cm 4 ). Thus it is very 
likely that C42H18 is completely dehydrogenated in the diffuse 
ISM and this could explain why it was not detected by Kokkin 
[eTaLl ( [2008] l. 

In the particular case of C20H4Q in the RR nebula, the situa- 
tion is less clear as the radiation field of this object is remarkably 
poor in UV photons, so that Fig. [10] cannot be used for this ob- 
ject. U sing the radiation field spectrum proposed by Mulas et al. 



2006JI, a density of «h = 10 cm with an arbitrary fraction 



of molecular hydrogen of 0.5, and a temperature of 100 K ( |Bu- 
jarra bal et al.|2005) , our model predicts coronene C2A^n to be 
completely dehydrogenated in RR. It is therefore very likely that 
corannulene C20H10 cannot survive in RR as well. As a conse- 
quence, our results strongly suggest that future search for PAH 
spectroscopic signatures should focus on much larger species, 
containing at least 50 or 60 carbon atoms. 

7.2. The carbon content in PAHs and carbon clusters 

The fact that we predict a large abundance of carbon clusters in 
the cavity and in a significant fraction of the PDR of NGC 7023 
NW raises the question of the evolution of carbon content in 
PAHs and carbon clusters. Carbon clusters have properties that 
differ significantly from those of PAHs, in terms of structure, 
stability and spectroscopy. Therefore, when a PAH is photodis- 
sociated until full dehydrogenation, it can be considered as a loss 
of carbon content in the PAH population and a gain for the car- 
bon cluster population. However, the latter population also may 
undergo evolution via photodissociation and carbon accretion. 
To get a first estimate, we used our model for computing the 
effective photodissociation rates (number of dissociation events 
per time unit) of planar C^" 1 ", and C°g + in the cavity and 
the NW PDR. We used the IR photon emission rates presented 



of graphite. 

We found that the largest cluster in this study, C^ 6 + , is de 
stroyed within a time longer than the age of the star, even at 
10" from the star. C"g + needs more than several 10 6 years to 
be significantly photodissociated in the NW PDR, and its de- 
struction rate rises steeply in the cavity until very fast photode- 
struction (~ 10 6 s _1 ) as multiphoton events become frequent for 
the corresponding radiation field. photodissociates much 

faster than larger species in the NW PDR. Its photodissociation 
rate increases less rapidly than for larger species when going 
closer to the star, and reaches the same value as for C^ 6 + close 
to the star. Here again, the evolution of the photodestruction 
rate is highly non-linear due to the predominance of multipho- 
ton events (threshold energies for fragmentation are ~15, 37, 54 
eV for C% + , and C^ + , respectively). 

These results draw a picture in which the cavity is populated 
by large carbon clusters while the PDR hosts a blend of large 
PAHs and relatively small carbon clusters. Carbon clusters are 
expected to explore various isomeric forms including non pla- 
nar (cage) ones (see for instance |Jones & Seifert|1997| for the 
various forms of C24). Heating by UV photons favours the ex- 
ploration of the various isomeric forms. Our results are therefore 
consistent with the formation of the cage fullerene C60 in the cav- 
ity following the destruction of PAHs as suggested by |Berne~&| 
|Tielens| ( |20fel l. 

In addition, this raises the question of the contribution of car- 
bon clusters to the AIBs, as well as to the physical evolution of 
PDRs, e.g. through the photoelectric effect. Further considera- 
tion on these questions calls for more laboratory studies on the 
photophysical and chemical properties of these species. 



7.3. Perspectives 

In this paper, we have focused on the interaction with photons 
and the reactivity with electrons and hydrogen. Including a 
larger set of chemical processes would be of interest for future 
studies. 

The reactivity of PAHs with C + has been studied experimen- 
tally by Cano sa~et al.| ( |1995) . The authors found the reaction to 
proceed at a rate close to the rate of Langevin with, in the case 
of anthracene, two possible channels of comparable branching 
ratio: charge exchange (c.e.) and carbon accretion. In Sect. |5.3| 
we have evaluated the timescale for c.e. to be ~ 8 years in the 
conditions of NGC 7023 NW. A similar timescale is expected for 
carbon accretion. These values are shorter than the timescale to 
reach hydrogenation steady state, which indicate that the cor- 
responding reactions have to be included in a PAH chemical 
model. Still, reactions of PAHs with C + have been studied only 
for small PAHs and significant variations with the size of PAHs 
are observed in the branching ratio between c.e. and carbon ac- 
cretion, calling for investigation of the reactivity of larger PAHs, 
including dehydrogenated PAHs, with C + . 

PAHs can also react with atomic oxygen but these reactions 
are also poorly known. |Le Page et al.| ( |1999a| l considered a reac- 
tion rate for PAH cations of fc+o ~ 10~ lu cm 3 s~' for the addition 
of an O-atom, regardless of the size of PAH cations. Assum- 
ing a density «h = 10 4 ctrT 3 and an O/H ratio of 3 x 10~ 4 , this 
leads to a reaction timescale of ~ 100 years. Therefore, for large 
PAHs, reaction with oxygen could play a significant role in PAH 
evolution. 
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Finally, it might be interesting to investigate the possibil- 
ity that PAHs form clusters and/or complexes with heavy atoms 
such as iron and silicon as proposed by several authors (Rapaci- 
|oli et al.|2006"l|Simon & Joblin|2009[|Joalland et al.|2009) i~ 

8. Conclusion 

In this paper, we presented a new model dedicated to the evo- 
lution of the hydrogenation and charge states of PAHs. This 
model was designed to compute the time evolution of the in- 
ternal energy of species, so that multi-photon dissociations, in 
the sense of the successive absorptions of several photons be- 
fore the species could completely cool down, are taken into ac- 
count for every dissociation process. We applied this model 
to four compact PAHs, namely coronene (C24H12), circum- 
coronene (C54H18), circumovalene (C66H20) and circumcircum- 
coronene (C96H24), in the conditions of NGC 7023, both in the 
NW PDR and the cavity, and of the diffuse ISM. 

The calculated spatial evolution of the PAH charge state is 
in good agreement with the observational constraints. The hy- 
drogenation state is more difficult to predict as it depends on 
a wealth of processes, many of which have not or poorly been 
studied in the laboratory or theoretically. Clear observational di- 
agnostics of the hydrogenation state of PAHs are also missing. 
Despite these difficulties, we showed that the general behaviour 
of PAH evolution is quite robust with molecular size and with the 
assumptions regarding the molecular processes. It leads to four 
major populations: superhydrogenated PAHs, normally hydro- 
genated PAHs, partially dehydrogenated PAHs and carbon clus- 
ters (i.e. totally dehydrogenated PAHs). 

Partially dehydrogenated PAHs are found to be minor 
species in terms of abundance regardless of the molecular size 
or the position in the PDR. Superhydrogenated PAHs dominate 
only for very large PAHs (Nc Z 96) in the conditions of the 
NGC 7023 NW PDR. The size limit below which normally hy- 
drogenated PAHs fully dehydrogenate to give birth to carbon 
clusters is found to lay between Nc ~ 50 and 70 at the emis- 
sion peak of AIBs, depending on the reaction rates of PAHs 
with H and H2. In the more diffuse medium of the cavity, where 
the fullerene C60 was recently observed, we predict all the four 
species studied here to be fully dehydrogenated. 

We evaluated that the electronic recombination of PAH 
cations with electrons and the reactivity of PAHs with H2 are 
the main contributors to the uncertainties on our results. 

In addition, we were able to rationalize the non detection of 
corannulene C20H10 in the Red Rectangle nebula and of hexa- 
peri-hexabenzocoronene 042^3 in the diffuse ISM. Therefore 
future search for the identification of individual PAHs should 
focus on larger species, bearing at least 50-60 carbon atoms. 

More generally, multi-photon events were found to dominate 
the evolution of interstellar PAHs, even for relatively low radi- 
ation fields as prevailing in the diffuse ISM. This results in a 
complex behaviour with variations of n H an d Go, that cannot be 
simply parametrised with the ratio Go/«h- 

We also determined PAH evolution timescales to be much 
longer than previously thought. In some cases, timescales com- 
parable to dynamical evolution timescales are reached, raising 
the question of the coupling between dynamics and chemical 
evolution of PAHs. This could impact both the physics and 
chemistry of PDRs, for example through photoelectric heating 
and H2 formation on PAHs. An other important consequence of 
these long timescales is the fact that reactions with species like O 
and C + may occur at similar or shorter timescales and therefore 
require more investigation. 



Finally, we showed that a large population of carbon clusters 
could exist in PDRs as a result of the dehydrogenation of PAHs. 
Our results therefore call for fundamental studies on carbon 
clusters in order to determine to which extent these species 
contribute to the mid-IR emission, and affect the physical and 
chemical evolution of PDRs. These studies are also necessary 
to unveil the link between PAHs and fullerenes. 
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Appendix A: Modelling the internal energy of PAHs 

In this work, we model the time evolution of both the abun- 
dances and internal energies of the studied species. The latter 
is described as a histogram in which each bin i represents the 
number per unit of volume of species bearing an internal energy 
between = iAE and = (z + 1)AE, where AE was set to 
0.25 eV. In this appendix, we detail how we relate the rates char- 
acterising the fluxes of species from one bin to another, which 
are driven by UV- visible photon absorption, IR photon emission, 
and molecule dissociation. 

Appendix A.1: Heating molecules 

Species are heated after the non-ionizing absorption of a UV- 
visible photon. This leads to a transition from one bin i to a 
higher bin i' (f > i), at a rate which depends on the flux of 
UV-visible photons T(hv), the absorption cross-section of the 
species <x abs (/iy), and the non-ionization yield 1 - Ti on (/iv). It also 
depends on the triangle shaped probability function Vi^ijiv). 
This function expresses the probability that the energy hv of the 
photon falls between the minimum (E?-i - Ef) and maximum 
(£V - Ei-i) differences between the boundaries of bins i and /'. 
Therefore, the rate of transitions from bin i to bin i' is given by: 

kl L= [ E " *V(*V> - Y ion (hv)) cr a Uhv)<F(hv)d(hv) 

(A.l) 

This expression is accurate, as long as one can neglect the varia- 
tions of o" a bs and Tic,, with the internal energy of the species, and 
the variations of abundance within the bin i. The latter condition 
requires AE to be smaller than hv for the less energetic absorbed 
photons that have a significant impact on the internal energy of 
the species. This is a much weaker constraint that the conditions 
for a proper description of the cooling of molecules. 
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Appendix A.2: Cooling molecules 

Species cool down through the emission of IR photons. This 
corresponds to the transition from one bin i' to a lower bin i (/' > 
z), at a rate that depends both on the IR emission rate and on the 
energy carried by the emitted photon which has to fall between 
the minimum and maximum differences between the boundaries 
of bins i and i' . We define the factor 6^ (E) that equals 1 when 
the latter condition is true, and otherwise. For the sake of 
simplicity, and to ensure the computations to be tractable, we 
neglect the dispersion of the IR energies of the emitted photons 
and compute Sjj>{E) using the mean IR energy as presented in 
Sect. 3.3.1 Therefore, the effective cooling rate from bin i' to 



bin i is given by: 



(E)k 1R (E)dE 



(A.2) 



This formulation requires to use bins thinner than the less ener- 
getic IR photon to ensure that each emission leads to an actual 
cooling. 

Appendix A.3: Molecule dissociation 

The dissociation rate from the bin i of a species is computed as 
the mean value of the dissociation rate of the species over the 
width AE of the bin: 



"diss 



JiAE 



(H-l)AE 



k diss (E) dE/AE 



(A.3) 



This approximate expression is better for smaller energy bins. 



Appendix A.4: Width of energy bins 

We empirically found that AE = 0.25 eV provides abundances 
accurate at ~ 5 % of the asymptotic results obtained when AE —> 
0. This value is therefore a good compromise between accuracy 
and computing time. 
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